Overcoming conventional modeling limitations using image- driven lattice-boltzmann method simulations for biophysical applications

The challenges involved in modeling biological systems are significant and push the boundaries of conventional modeling. This is because biological systems are distinctly complex, and their emergent properties are results of the interplay of numerous components/processes. Unfortunately, conventional modeling approaches are often limited by their inability to capture all these complexities. By using in vivo data derived from biomedical imaging, image-based modeling is able to overcome this limitation. In this work, a combination of imaging data with the Lattice-Boltzmann Method for computational fluid dynamics (CFD) is applied to tissue engineering and thrombogenesis. Using this approach, some of the unanswered questions in both application areas are resolved. In the first application, numerical differences between two types of boundary conditions: “wall boundary condition” (WBC) and “periodic boundary condition” (PBC), which are commonly utilized for approximating shear stresses in tissue engineering scaffold simulations is investigated. Surface stresses in 3D scaffold reconstructions, obtained from high resolution microcomputed tomography images are calculated for both boundary condition types and compared with the actual whole scaffold values via image-based CFD simulations. It is found that, both boundary conditions follow the same spatial surface stress patterns as the whole scaffold simulations. However, they under-predict the absolute stress values approximately by a factor of two. Moreover, it is found that the error grows with higher scaffold porosity. Additionally, it is found that the PBC always resulted in a lower error than the WBC. In a second tissue engineering study, the dependence of culture time on the distribution and magnitude of fluid shear in tissue scaffolds cultured under flow perfusion is investigated. In the study, constructs are destructively evaluated with assays for cellularity and calcium deposition, imaged using µCT and reconstructed for CFD simulations. It is found that both the shear stress distributions within scaffolds consistently increase with culture time and correlate with increasing levels of mineralized tissues within the scaffold constructs as seen in calcium deposition data and µCT reconstructions. In the thrombogenesis application, detailed analysis of time lapse microscopy images showing yielding of thrombi in live mouse microvasculature is performed. Using these images, image-based CFD modeling is performed to calculate the fluid-induced shear stresses imposed on the thrombi’s surfaces by the surrounding blood flow. From the results, estimates of the yield stress (A critical parameter for quantifying the extent to which thrombi material can resist deformation and breakage) are obtained for different blood vessels. Further, it is shown that the yielding observed in thrombi occurs mostly in the outer shell region while the inner core remains intact. This suggests that the core material is different from the shell. To that end, we propose an alternative mechanism of thrombogenesis which could help explain this difference. Overall, the findings from this work reveal that image-based modeling is a versatile approach which can be applied to different biomedical application areas while overcoming the difficulties associated with conventional modeling

Global, regional, and national incidence and mortality for HIV, tuberculosis, and malaria during 1990–2013: a systematic analysis for the Global Burden of Disease Study 2013

BACKGROUND: The Millennium Declaration in 2000 brought special global attention to HIV, tuberculosis, and malaria through the formulation of Millennium Development Goal (MDG) 6. The Global Burden of Disease 2013 study provides a consistent and comprehensive approach to disease estimation for between 1990 and 2013, and an opportunity to assess whether accelerated progress has occured since the Millennium Declaration. METHODS: To estimate incidence and mortality for HIV, we used the UNAIDS Spectrum model appropriately modified based on a systematic review of available studies of mortality with and without antiretroviral therapy (ART). For concentrated epidemics, we calibrated Spectrum models to fit vital registration data corrected for misclassification of HIV deaths. In generalised epidemics, we minimised a loss function to select epidemic curves most consistent with prevalence data and demographic data for all-cause mortality. We analysed counterfactual scenarios for HIV to assess years of life saved through prevention of mother-to-child transmission (PMTCT) and ART. For tuberculosis, we analysed vital registration and verbal autopsy data to estimate mortality using cause of death ensemble modelling. We analysed data for corrected case-notifications, expert opinions on the case-detection rate, prevalence surveys, and estimated cause-specific mortality using Bayesian meta-regression to generate consistent trends in all parameters. We analysed malaria mortality and incidence using an updated cause of death database, a systematic analysis of verbal autopsy validation studies for malaria, and recent studies (2010-13) of incidence, drug resistance, and coverage of insecticide-treated bednets. FINDINGS: Globally in 2013, there were 1·8 million new HIV infections (95% uncertainty interval 1·7 million to 2·1 million), 29·2 million prevalent HIV cases (28·1 to 31·7), and 1·3 million HIV deaths (1·3 to 1·5). At the peak of the epidemic in 2005, HIV caused 1·7 million deaths (1·6 million to 1·9 million). Concentrated epidemics in Latin America and eastern Europe are substantially smaller than previously estimated. Through interventions including PMTCT and ART, 19·1 million life-years (16·6 million to 21·5 million) have been saved, 70·3% (65·4 to 76·1) in developing countries. From 2000 to 2011, the ratio of development assistance for health for HIV to years of life saved through intervention was US$4498 in developing countries. Including in HIV-positive individuals, all-form tuberculosis incidence was 7·5 million (7·4 million to 7·7 million), prevalence was 11·9 million (11·6 million to 12·2 million), and number of deaths was 1·4 million (1·3 million to 1·5 million) in 2013. In the same year and in only individuals who were HIV-negative, all-form tuberculosis incidence was 7·1 million (6·9 million to 7·3 million), prevalence was 11·2 million (10·8 million to 11·6 million), and number of deaths was 1·3 million (1·2 million to 1·4 million). Annualised rates of change (ARC) for incidence, prevalence, and death became negative after 2000. Tuberculosis in HIV-negative individuals disproportionately occurs in men and boys (versus women and girls); 64·0% of cases (63·6 to 64·3) and 64·7% of deaths (60·8 to 70·3). Globally, malaria cases and deaths grew rapidly from 1990 reaching a peak of 232 million cases (143 million to 387 million) in 2003 and 1·2 million deaths (1·1 million to 1·4 million) in 2004. Since 2004, child deaths from malaria in sub-Saharan Africa have decreased by 31·5% (15·7 to 44·1). Outside of Africa, malaria mortality has been steadily decreasing since 1990. INTERPRETATION: Our estimates of the number of people living with HIV are 18·7% smaller than UNAIDS's estimates in 2012. The number of people living with malaria is larger than estimated by WHO. The number of people living with HIV, tuberculosis, or malaria have all decreased since 2000. At the global level, upward trends for malaria and HIV deaths have been reversed and declines in tuberculosis deaths have accelerated. 101 countries (74 of which are developing) still have increasing HIV incidence. Substantial progress since the Millennium Declaration is an encouraging sign of the effect of global action. FUNDING: Bill & Melinda Gates Foundation

Time-Dependent Shear Stress Distributions during Extended Flow Perfusion Culture of Bone Tissue Engineered Constructs

Flow perfusion bioreactors have been extensively investigated as a promising culture method for bone tissue engineering, due to improved nutrient delivery and shear force-mediated osteoblastic differentiation. However, a major drawback impeding the transition to clinically-relevant tissue regeneration is the inability to non-destructively monitor constructs during culture. To alleviate this shortcoming, we investigated the distribution of fluid shear forces in scaffolds cultured in flow perfusion bioreactors using computational fluid dynamic techniques, analyzed the effects of scaffold architecture on the shear forces and monitored tissue mineralization throughout the culture period using microcomputed tomography. For this study, we dynamically seeded one million adult rat mesenchymal stem cells (MSCs) on 85% porous poly(l-lactic acid) (PLLA) polymeric spunbonded scaffolds. After taking intermittent samples over 16 days, the constructs were imaged and reconstructed using microcomputed tomography. Fluid dynamic simulations were performed using a custom in-house lattice Boltzmann program. By taking samples at different time points during culture, we are able to monitor the mineralization and resulting changes in flow-induced shear distributions in the porous scaffolds as the constructs mature into bone tissue engineered constructs, which has not been investigated previously in the literature. From the work conducted in this study, we proved that the average shear stress per construct consistently increases as a function of culture time, resulting in an increase at Day 16 of 113%

Scaffolds with a High Surface Area-to-Volume Ratio and Cultured Under Fast Flow Perfusion Result in Optimal O2 Delivery to the Cells in Artificial Bone Tissues

Tissue engineering has the potential for repairing large bone defects, which impose a heavy financial burden on the public health. However, difficulties with O2 delivery to the cells residing in the interior of tissue engineering scaffolds make it challenging to grow artificial tissues of clinically-relevant sizes. This study uses image-based simulation in order to provide insight into how to better optimize the scaffold manufacturing parameters, and the culturing conditions, in order to resolve the O2 bottleneck. To do this, high resolution 3D X-ray images of two common scaffold types (salt leached foam and non-woven fiber mesh) are fed into Lattice Boltzmann Method fluid dynamics and reactive Lagrangian Scalar Tracking mass transfer solvers. The obtained findings indicate that the scaffolds should have maximal surface area-to-solid volume ratios for higher chances of the molecular collisions with the cells. Furthermore, the cell culture media should be flown through the scaffold pores as fast as practically possible (without detaching or killing the cells). Finally, we have provided a parametric sweep that maps how the molecular transport within the scaffolds is affected by variations in rates of O2 consumption by the cells. Ultimately, the results of this study are expected to benefit the computer-assisted design of tissue engineering scaffolds and culturing experiments